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Quantum dynamics (e.g., the Schrodinger equation) and classical dynamics (e.g., Hamilton equations) can 
both be formulated in equal geometric terms: a Poisson bracket defined on a manifold. The difference between 
both worlds is due to the presence of extra structure in the quantum case, that leads to the appearance of 
the probabilistic nature of the measurements and the indetermination and superposition principles. In this 
paper we first show that the quantum-classical dynamics prescribed by the Ehrenfest equations can also be 
formulated within this general framework, what has been used in the literature to construct propagation 
schemes for Ehrenfest dynamics. Then, the existence of a well defined Poisson bracket allows to arrive to a 
Liouville equation for a statistical ensemble of Ehrenfest systems. The study of a generic toy model shows 
that the evolution produced by Ehrenfest dynamics is ergodic and therefore the only constants of motion are 
functions of the Hamiltonian. The emergence of the canonical ensemble characterized by the Boltzmann's 
distribution follows after an appropriate application of the principle of equal a priori probabilities to this case. 
This work provides the basis for extending stochastic methods to Ehrenfest dynamics. 



I. INTRODUCTION 

There is little doubt that the solution of Schrodinger 
equation for a combined system of electrons and nuclei 
enables us to predict most of the chemistry and molecu- 
lar physics that surrounds us, including bio-physical pro- 
cesses of great complexity. Yet there is also little doubt 
that this task is not possible in general, and approxima- 
tions need to be made; one of the most important and 
successful being the classical approximation for a num- 
ber of the particles. Mixed quantum-classical dynamical 
(MQCD) models are therefore necessary and widely used. 

The so-called "Ehrenfest equations" (EE) result from 
a straightforward application of the classical limit to a 
portion of the particles of a full quantum system, and 
constitute the most evident MQCD model, as well as a 
first step in the intricatcd problem of mixing quantum 
and classical dynamics. For example, it can be noted 
that much of the field of molecular dynamics (MD) - 
whether ab initio or not - is based on equations that can 
be obtained from further approximations to the Ehrenfest 
model: for instance, written in the adiabatic basis, the 
Ehrenfest dynamics collapse into Born Oppcnheimer MD 
if we assume the non-adiabatic couplings to be negligible. 

It is not the purpose of this work, however, to dwell into 
the justification or validity conditions of the EE (see for 
example Refs. [lH^ for rigurous analyses). Nor is it to in- 
vestigate the unsettled problem of which is the best man- 
ner of mixing quantum and classical degrees of freedom. 
The Ehrenfest model has a proven application niche, and 
for this reason we are interested in investigating some of 
its theoretical foundations, in a manner and with an aim 
that we describe in the following^ 

Classical mechanics (CM) can be formulated in several 
mathematical frameworks, each corresponding to a differ- 



ent level of abstraction (Newton's equations, the Hamil- 
tonian formalism, the Poisson brackets, etc.). Perhaps 
its more abstract and general formulation is geometri- 
cal, in terms of Poisson manifolds. Similarly, quantum 
mechanics (QM) can be formulated in different ways, 
some of which resemble its classical counterpart. For 
example, the observables (self-adjoint linear operators) 
are endowed with a Poisson algebra almost equal to 
the one that characterizes the dynamical variables in 
CM. Moreover, Schrodinger equation can be recast into 
Hamiltonian equations form 7 - by transforming the com- 
plex Hilbert space into a real one of double dimension; 
the observables are also transformed into dynamical func- 
tions in this new phase space, in analogy with the classi- 
cal case. Finally, a Poisson bracket formulation has also 
been established for QM, which permits to classify both 
the classical and the quantum dynamics under the same 
heading. 

This variety of formulations does not emerge from aca- 
demic caprice; the successive abstractions simplify fur- 
ther developments of the theory, such as the step from 
microscopic dynamics to statistical dynamics: the deriva- 
tion of Liouville's equation (or von Neumann's equation 
in the quantum case), at the heart of statistical dynam- 
ics, is based on the properties of the Poisson algebra. 

The issue regarding what is the correct equilibrium 
distribution of a mixed quantum-classical system is seen 
as a very relevant one 4^— An attempt to its derivation 
can be found in Ref. Il4l . where they arrive to the same 
distribution that we will advocate below, although it is 
done using the Nose-Hoover technique , 16 ' 17 which is only 
a mathematical scheme to produce the equilibrium, and 
which the very authors of Ref. agree that it is not clear 
how to apply to a mixed quantum-classical dynamics. In 
Ref. [H they provide some analytical results about this 
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distribution, but not in the case in which the system of 
interest is described by the EE; instead (and as in Refs.H 
and [13), they assume that the system is fully quantum, 
and that it is coupled to an infinite classical bath via an 
Ehrenfest-like interaction. 

It is therefore necessary to base Ehrcnfcst dynam- 
ics - or any other MQCD model - on firm theoretical 
ground. In particular, we are interested in establishing 
a clear path to statistical mechanics for Ehrenfest sys- 
tems, which in our opinion should be done by first embed- 
ding this dynamics into the same theoretical framework 
used in the pure classical or quantum cases (i.e., Poisson 
brackets, symplectic forms, etc). Then the study of their 
statistics will follow the usual steps for purely classical 
or purely quantum ensembles. In this respect, it should 
be noted that other approaches to MQCD (not based on 
Ehrenfest equations) exist ^r— and to their correspond- 
ing statistics) 21 ' 22 however it has been found that the for- 
mulation of well defined quantum-classical brackets (i.e., 
satisfying the Jacobi identity and Leibniz derivation rule) 
is a difficult issued - — On the contrary, as shown below, 
and as a result of the fact that the evolution of a quantum 
system can be formulated in terms of Hamilton equations 
similar to those of a classical system^21 we will not have 
difficulties in a rigorous formulation of the Ehrenfest dy- 
namics in terms of Poisson brackets. 

The roadmap of this project is the following: In Sec- 
tion|TT]we recall the definition of the Ehrenfest model. In 
Section irm we summarize quickly the formulation of CM 
in terms of Poisson brackets. Then we summarize, in 
Section IIV1 the analogous description of QM in terms 
of geometrical objects which can be found, for exam- 
ple, in Refs. 0, I28U351 These works demonstrate how 
Schrodinger equation can be written as a set of (appar- 
ently classical) Hamiltonian equations. Also, by using 
a suitable definition of observables as functions on the 
(real) set of physical states, Schrodinger equation can be 
written in terms of a canonical Poisson bracket. QM 
appears in this way as "classical", although the exis- 
tence of extra algebraic structure encodes the probabilis- 
tic interpretation of measurements, and the superposi- 
tion and indetermination principles 2 ^— . For the case of 
Ehrenfest dynamics, we can now combine it (Section [Vj 
with CM, whose Poisson bracket formulation is very well 
known. The procedure is very simple: we consider as 
a global phase space the Cartesian product of the elec- 
tronic and the nuclear phase spaces and define a global 
Poisson bracket simply as the sum of the classical and the 
quantum ones. It is completely straightforward to prove 
that such a Poisson bracket is well defined and that it 
provides the correct dynamical equations (EE). 

From a formal point of view, the resulting dynamics 
is more similar to a classical system than to a quantum 
one, although when considering a pure quantum system 
the dynamics is the usual Schrodinger one. This is not 
surprising since the coupling of the classical and quantum 
systems makes the total system nonlinear in its evolution, 
and this is one of the most remarkable differences between 



a classical and a quantum system in any formulation. 

Once the dynamical description as a Poisson system 
is at our disposal, it is a very simple task (Section IVI|) 
to construct the corresponding statistical description fol- 
lowing the lines of Refs. HH and[37]. The formal similarity 
with the CM case ensures the correctness of the proce- 
dure and allows us to derive a Liouville equation. The 
choice of the equilibrium distribution is based, as usual, 
in the principle of equal a priori probablities implicitely 
used by Gibbs and clearly formulated by Tolman^. 

Finally, in the conclusions (Section IVII[) . we propose 
to extend, within our formalism, stochastic methods to 
Ehrenfest dynamics. 



II. THE EHRENFEST MODEL 

The Ehrenfest Equations have the following general 
form: 



Rj(t) = m) 
pj(t) = 



OH 
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(R(t),P(t),t)m)), (2-1) 

(R(t),P(t),t)m)), (2-2) 
(2.3) 



ih-m)) = H(R(t),p(t),t)m)) 



where (R, P) denote collectively the set of canonical po- 
sition and momenta coordinates of a set of classical parti- 
cles, whereas ip is the wavefunction of the quantum part 
of the system. See Ref [3^ for the issue of the Hellmann- 
Feynman theorem in this context - i.e., whether one 
should take the derivative inside or outside the expec- 
tation value. 

This and other MQCD models appear in different con- 
texts; in many situations, the division into quantum and 
classical particles is made after the electrons have been 
integrated out, and is used to quantize a few of the nu- 
clear degrees of freedom. However, a very obvious case 
to use EE is when we want to treat electrons quantum 
mechanically, and nuclei classically. In order to fix ideas, 
let us use this case as an example: the Hamiltonian for 
the full quantum system is: 
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(2.4) 



where all sums must be understood as running over the 
whole natural set for each index. Mj is the mass of the 
J-th nucleus in units of the electron mass, and Zj is the 
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charge of the J-th nucleus in units of (minus) the elec- 
tron charge. Also note that we have defined the nuclei- 
electrons potential V n - e (r,R) and the electronic Hamil- 
tonian H e (r,R) operators. 

The EE may then be reached in the following way>i~— 
first, the full wave function is split into a product of nu- 
clear and electronic wave functions, which leads to the 
time dependent self-consistent field model, in which the 
two subsystems arc quantum and coupled. Afterwards, 
a classical limit procedure is applied to the nuclear sub- 
system, and the EE emerge naturally. In terms of the 
nuclei positions Rj and of the element of the Hilbert 
space G % which encodes the state of the electrons 
of the system, the Ehrenfest dynamics is then given by 



MjRj 



-(iP\VjH e (r,R)\iP), 



ih-\il>)=H e {r,R)\i)). 



(2.5) 
(2.6) 



These equations can be given a Hamiltonian-type de- 
scription by introducing a Hamiltonian function of the 
form: 



P 



H(R,P) = ^ + WHeF, RM). 



(2.7) 



Then, by fixing a relation of the form Pj = MRj, we 
obtain a structure similar to Hamilton equations: 



Rj 



Pj_ 

My 



ih- 



dr 



Pj = -(TP\VjH e (r,R)\iP) ) 
= H e (f,R)\ip). 



(2.8) 

(2.9) 
(2.10) 



But in spite of the formal similarities, these equations do 
not correspond yet to Hamilton equations, since they lack 
a global phase space formulation (encompassing both the 
nuclear and the electronic degrees of freedom) and a Pois- 
son bracket. 



III. CLASSICAL MECHANICS IN TERMS OF POISSON 
BRACKETS 

Let us begin by recalling very quickly the Hamiltonian 
formulation of classical dynamics. We address the inter- 
ested reader to a classical text, such as Ref. H(| for a 
more detailed presentation. 

Let us consider a classical system with phase space 
Mc, which, for the sake of simplicity, we can identify 
with K 2 ™, where n is the number of degrees of freedom 
(strictly speaking, Mc is a general 2n-dimcnsional mani- 
fold, homeomorphic to M. 2n only locally). The 2n dimen- 
sions correspond with the n position coordinates that 
specify the configuration of the system, and the n cor- 
responding momenta (mathematically, however, the di- 
vision into "position" and "momenta" coordinates is a 
consequence of Darboux theorem - see below). 



In what regards the observablcs, classical mechanics 
uses the set of diffcrcntiablc functions 



f-M G 



(3.1) 



assigning the result of the measurement to every point 
in Mc- On this set of functions C°°(Mc) we introduce 
an operation, known as Poisson bracket, which allows 
us to study the effect of symmetry transformations and 
also the dynamical evolution. The precise definition is as 
follows: 

Definition 1 A Poisson bracket, {■,■}, is a bilinear 
operation 

{-, •} : C°°{M C ) x C°°{M C ) -> C°°{M C ), (3.2) 
which: 

• It is antisymmetric, 

{/,<?} = -{<?,/}, V/, 5 er(M c ), 

• It satisfies the Jacobi identity, i.e., \ff,g,h £ 
C°°{M C ): 

{/, {<?, h}} + {h, {f,g}} + {g, {h, /}} = 0. 

• It satisfies the Leibniz rule ,i.e., V/, g, h £ 
C°°(M C ): 

{f,gh} = {f,g}h + g{f,h}. 

A Poisson bracket allows us also to introduce the con- 
cept of Hamiltonian vector field: 

Definition 2 Given a function f S C°°(Mc) and a 
Poisson bracket {•,•}, a vector field, Xf, is said to be 
its Hamiltonian vector field if 

X f (g) = {f,g}, VgeC°°(M G ). 

Definition 3 We shall call a Hamiltonian system to 

a triple (Mc, {■,■}, H), where {•,•} is a Poisson bracket 
on Mc, and dynamics is introduced via the function H £ 
C°°(Mc), that we call the Hamiltonian. 

If the Poisson bracket is non degenerate (it has no 
Casimir functions), a theorem due to Darboux ensures 
that there exists a set of coordinates (R, P) for which 
the bracket has the "standard" form (at least locally, in 
a neighbourhood of every point): 



dR k dP k 



f-' dP k dR k 

fe=i K 



(3-3) 



These coordinates are called Darboux coordinates. They 
are specially useful when studying the dynamics and in- 
variant measures of a Hamiltonian systems!*. In the rest 
of the paper we shall be working with this kind of coor- 
dinates. 
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Then, given / 6 C°°(M. 2n ), we can write the corre- 
sponding Hamiltonian vector field Xf as 



X 



f 



df(R,P) d df(R,P) d 



^ dP k dR k 

k=l 



dR k dP k ' 



(3.4) 



Now one can consider two different formulations of the 
dynamics: 

• One which defines the corresponding Hamiltonian 
vector field Xh obtained as above 

X H (g) = {H,g}, VgeC°°(M c ). 

The integral curves of the vector field Xh define 
the solution of the dynamics. 

• An analogous formulation can be given in terms 
of the observablcs. If we consider now the set of 
functions of the system, i.e., the set of classical ob- 
scrvables, the dynamics is written as the Poisson 
bracket of the Hamiltonian function H with any 
other function of the system, i.e., 

f| = {#,/}, V/eC°°(A/ c ). (3.5) 

Both approaches are equivalent: the differential equa- 
tions that determine the integral curves of the Hamilto- 
nian vector field are given by Eq. 13.51 for the functions 
"position" and "momenta" of each particle (i.e., R and 
P) . These equations are nothing else but Hamilton equa- 
tions: 

dH 



W = 



P, 



dH 

dRi 



(3.6) 



IV. SUMMARY OF GEOMETRIC QUANTUM 
MECHANICS 

The aim of this section is to provide a description of a 
quantum mechanical system by using the geometric tools 
which are used to describe classical mechanical systems. 
It is just a very quick summary of the framework which 
has been developed in the last 30 years and which can 
be found in Refs. 0, [2814351 and references therein. For 
the sake of simplicity, we shall focus only on the finite 
dimensional case. The Hilbert space T-L becomes then 
isomorphic to C n for n a natural number. 



A. The states 

Consider a basis {IV'fe}} f° r H- Each state \ip) £ H. can 
be written in that basis with complex components (or 
coordinates, in more differential geometric terms) {z k }: 



H3 M = 5>kllM- 



We can just take the vector space inherent to the Hilbert 
space, and turn it into a real vector space Mq, by split- 
ting each coordinate in its real and imaginary part: 



%3 z k = q k + ip k h-> (cj k ,p k ) G 



1)2/ 1 



M, 



Q- 



We will use real coordinates (q k ,p k ), k = 1, . . . ,n to rep- 
resent the points of T-L when thought as real manifold 
elements. From this real point of view the similarities 
between the quantum dynamics and the classical one 
described in the previous section will be more evident. 
Sometimes it will be useful to maintain the complex no- 
tation tjj or z k for the elements of the Hilbert space. 

Another important aspect of the Hilbert space descrip- 
tion of quantum mechanics is the study of the global 
phase of the state. It is a well known fact that physical 
states are independent from the global phase of the el- 
ement of the Hilbert space that we choose to represent 
them. In the formulation as a real vector space, we can 
represent the multiplication by a phase on the manifold 
Mq as a transformation whose infinitesimal generator is 
written as: 



z k 



8 
dz k 



d 

dz k 



E 



_d_ 

dpk 







Pk- 



(4.1) 

where the derivatives with respect to the complex and 
real variables are related by S— = h(ir i?r-) and 

J az k 2 v c% dp k > 



Szfc 2 V dq k 



Hi 



J dp k 



). The meaning of this vector field is 
simple to understand if we realize that a phase change 
modifies the angle of the complex number representing 
the state, when considered in polar form (i.e., in polar 



coordinates {r k ,9 k } 



fe=l, 



with z k 



rke 



Eq. (IXTD 



becomes T = ^2 k de k ). Then, from a geometrical point 
of view we can use Eq. (|4.1j) in two ways: 

• Computing its integral curves, which arc the differ- 
ent states which are obtained from an initial one 
by a global phase multiplication. 

• Acting with the vector field on functions of Mq 
(which will represent our observablcs) providing us 
with the effect of the global phase transformation 
on the observables. 

One final point is to consider the limitation in the 
norm. The sphere of states with norm equal to one in 
C™ corresponds in the real-vector-space description to 
the (2Nq — 1) -dimensional sphere: 



Sq = G Mq\ Y,{qt+pl) = 1 j 



(4.2) 



It is immediate that the vector field (|4.1|) is tangent to 
Sq, since the phase change preserves the norm of the 
state. 
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B. The Hermitian structure 

The Hilbert space structure of a quantum system is 
encoded in a Hermitian structure or scalar product (•, •}. 
The latter is specified by the choice of an orthonormal 
basis {\<f>k)}- 

Taking as coordinates z k = q k +ip k , the components of 
a vector in this basis, as we did in the previous paragraph, 
we can define a Poisson bracket in TL, by 

^\dz k dz k dzkdzkj 

or in real coordinates 

f , ! ( 9f dg df dg\ 

^ \ dp k dq k dq k dp k J 

that corresponds to the standar Poisson bracket in clas- 
sical mechanics. This justifies the choice of the notation 
for the real coordinates wich become Darboux coordi- 
nates for the Poisson bracket. 

For later purposes it will be useful to introduce a sym- 
metric bracket by 

^\dz k dz k dzkdzkj 

=E(i i i^+# L # L y (^) 

V d 9k OQk op k dp k J 

Observe that in order to define the Poisson and the sym- 
metric brackets we made use of an orthonormal basis in 
the Hilbert space. We may think that both objects are 
not canonical as they could depend of the choice of the 
basis. In the following proposition we shall show that 
this is not the case. 

Proposition 1 The brackets { , } and { , } + defined 
above depend only on the scalar product and not on the 
orthonormal basis used to construct them. 

Proof The coordinates in two different orthonormal basis, 
for a given scalar product, {IV'fc)} an d HV'fe)} are related 
by a unitary transformation, i.e., if 

k k 

we have z] = J2 k U?z k with J2 k = S ol- 

Therefore 

W~ = E u jJT7 and 7F~ = E jj=7 • 
dz k ^ ■> dz\ dz k ^-r 1 3 oz'. 

From which one immediately gets 

df_dg_ = df dg 
^ dz k dz k Z-s g z > Qz> 

k 3 3 3 



and both the Poisson and symmetric brackets are inde- 
pendent from the choice of orthonormal basis, which con- 
cludes the proof. □ 
An important property of the brackets defined above 
is that they are preserved by the vector field V in Eq. 
(|4.ip in the sense that 

r{/,. 9 } = {r/,.g} + {/,r 5 } 

and 

r{/,. g } + = {r/,. g } + + {/,r.g} + . 

A fact that can be proved with a simple computation. 

This property is important for us because it implies 
that the symmetric or antisymmetric bracket of two func- 
tions killed by Y (i.e., constant under phase change) will 
be killed by Y too: 

Yf = Yg = => g} = Y{f, g}+ = 0. 

Therefore, it makes sense to restrict ourselves to this class 
of functions. 

C. The observables 

Now that we have been able to formulate quantum 
mechanics on a real state space, we proceed to discuss 
how to represent the physical observables on this new 
setting. Instead of considering the observables as linear 
operators (plus the usual requirements, sclf-adjointncss, 
boundedness, etc.) on the Hilbert space W, we shall be 
representing them as functions defined on the real man- 
ifold Mq. The reason for that is to resemble, as much 
as possible, the classical mechanical approach. But we 
can not forget the linearity of the operators, and thus 
the functions must be chosen in a very particular way. 
The usual choice inspired in Ehrenfest's theorem selects 
the functions introduced in the following: 

Definition 4 To any operator A e Lin(H) we associate 
the quadratic function 

f A (4,) = ^\A^. (4.7) 

We shall denote the set of such functions as J- s {Mq). 

Notice that this definition of observable is different 
from the analogous one in the classical case. In the 
classical framework a state, represented by a point in 
Mq, provides a well defined result for any observable 
/ : Mq — > M. In this geometric quantum mechanics 
(GQM in the following) framework, on the other hand, 
the value at a given state of f A € J- s (Mq) provides just 
the average value of the corresponding operator in that 
state. Besides, contrarily to the classical case, not all the 
C°° functions on Mq are regarded as observables, but 
only those of the form of (|4.7[) . 
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Once the definition has been stated, we must verify 
that it is a consistent one, in the sense that the brackets 
defined in H act correctly in this space. This is contained 
in the following 

Lemma 1 The set of functions of the form (|4.7|) is 
closed under the brackets defined in section B. 

Proof Direct computation: just take the expression of the 
functions, compute the derivatives, obtain the bracket 
and verify that they correspond to the expressions of the 
functions associated to the (anti)commutator of the cor- 
responding operators. □ 
Within the usual approach of quantum mechanics, 
there are three algebraic structures on the set of oper- 
ators, which turn out to be meaningful and important 
for the physical description: 

• The associative product of two operators: 

A,B e Lin(H) A.B G Lin(H). 

It is important to notice, though, that this opera- 
tion is not internal in the set of Hcrmitian opera- 
tors (i.e., those associated to physical magnitudes), 
since the product of two Hermitian operators is not 
Hermitian, in general. 

• The anticommutator of two operators: 

A, Be Lm(H) -> [A,B]+ G Lin(H). 

• The commutator of two operators: 

A, Be Lin(H) -> i[A,B] e Lm(H). 

Notice that these last two operations are internal in the 
space of Hermitian operators. How do we translate these 
operations into the GQM scheme? 

To answer this question we recall the brackets de- 
fined in the previous subsection B. If we take G 
F.(Mq), 

• the anticommutator of two operators, becomes the 
symmetric bracket or Jordan product of the func- 
tions (see Ref. IZll) : 

[A,B}+ f [A , B]+ = {fA,f B }+ G T S (M Q ). 

• The commutator of two operators, transkates into 
the Poisson bracket of the functions: 

i[A,B] -> fi[A,B] = VaJb} G T s (Mq). 

Conclusion: The operators and their algebraic struc- 
tures are encoded in the set of functions with the brackets 
associated to the Hermitian product. 

Thus we see why it makes sense to consider only that 
specific type of functions: it is a choice which guaran- 
tees to maintain all the algebraic structures which are 
required in the quantum description. 



Another important property of the set of operators of 
quantum mechanics is the corresponding spectral theory. 
In any quantum system, it is of the utmost importance 
to be able to find eigenvalues and eigenvectors. We can 
summarize the relation between these objects and those 
in the new GQM scheme in the following result: 

Lemma 2 Let f a be the function associated to the ob- 
servable A. Then, if we consider ip ^ (The case ip = 
is not meaningful from the physical point of view), 

• the eigenvectors of the operator A coincide with the 
critical points of the function f a, i-e., 

df A (ip) = =>■ ip is an eigenvector of A. 

• The eigenvalue of A at the eigenvector ip is the 
value that the function f A takes at the critical point 
i>. 

Proof Trivial computation. Recall Ritz theorem from 
any quantum mechanics textbook, as for instance Ref. |42|. 
□ 

In the previous subsection we saw that the vector field 
defined in Eq. (|4.1|) ^ preserves the Poisson bracket. In- 
deed, it is possible to prove that T corresponds to a 
Hamiltonian vector field whose Hamiltonian function is 
precisely the norm of the state: 

Lemma 3 Let Mq = R 2 " . The vector field defining the 
multiplication of the states by a phase is the Hamiltonian 
vector field of the function corresponding to the identity 
operator: 

h = \wi>) — >x A = r. (4.8) 

Proof The function f\ is easily computed: 

k 

The Hamiltonian vector field corresponding to this func- 
tion is, easily too, computed: 




□ 



D. The dynamics 

As in the classical case, the dynamics can be imple- 
mented in different forms, always in a way which is com- 
patible with the geometric structures introduced so far: 

• In the Schrodinger picture, the dynamics is de- 
scribed as the integral curves of the vector field 
Xf H , where H is the Hamiltonian operator of the 
system: 

X fH =h- l {f H ,-}. (4.9) 



• In the Heiscnbcrg picture, the dynamics is intro- 
duced by transferring the Heiscnbcrg equation into 
the language of functions: 



fA = h- 1 {f H ,f A }. 



(4.10) 



The most interesting example is probably the expres- 
sion of the von Neumann equation: 



f P = fr l {f B ,f p }, 



(4.11) 



where f p is the function associated to the density matrix. 

Conclusion: It is possible to describe quantum dy- 
namics as a flow of a vector field on the manifold Mq or 
on the set of quadratic functions of the manifold. Such a 
vector field is a Hamiltonian vector field with respect to 
the Poisson bracket encoded in the Hermitian structure 
of the system. 

Example 1 Let us consider again the simplest quantum 
situation defined on C 2 . As a real manifold, Mq ~ R 4 . 
Consider then a Hamiltonian H : C 2 — > C 2 which is 
usually written as a matrix: 



H 



H21 H22 



If we consider it as a matrix on the real vector space Mq, 
it reads: 



Hjs = 



/ H qiqi H qipi H qiq2 H qip2 

H pigi H plPl H piq2 H plP2 

Hqiq\ -"92P1 Hqiq-i ^qivi 

^Hp 2 q 1 H p2pi H p2 q 2 H p2p2 / 



where the matrix above is symmetric because H is Her- 
mitian, since we have: 



H, 



qiqi ~ "IX — ±J -pipn 



H n 



Hq 



= H p 
= H rj 



H qiQ2 = Re(i?i2j = H PlP2 , H qiP2 = — lm{Hi2) = —H Piq2 , 
Hq 2Ql = Rc(i?2i) = H P2Pl , H q2Pl = — Im(i?2x) = —H P2qi . 

The function fjj in J- s (Mq) becomes thus: 

fH = l;ilA H KipR, where ifa = (qi,pi, ^2,^2)* 
and then, the Hamiltonian vector field turns out to be: 

fdf H d df H d , df H d df H d 



x H = nr 



\ dpi dqi dqi dpi dp2 dq 2 dq 2 dp-. 



and its integral curves are precisely the expression of 
Schrddinger equation when we write it back in complex 
terms: 

?x = fi~ (H piqi qi + H Pl q 2 q 2 + H PlPl pi + H PlP2 p 2 ), 
Pi = -ft~ 1 (H qi q 1 qi + H qiq2 q 2 + H qiPl pi + H qiP2 p 2 ), 
92 = h~ (H P2 q 2 q 2 + H P2qi qi + H P2P2 p 2 + H P2Pl pi), 

P2 = -h~ 1 (H q2 q 2 q 2 + Hq 2qi qi + H q2p2 p 2 + H q2pi pi). 



We can write these equations as: 
with 



(4.12) 



J = 



ro -1 
1000 
000-1 

\0 1 



or, equivalently, 



where 



\ip) = -h^iH^j) 



qi + ipi 
q_2 + ip2 



(4.13) 



(4.14) 



is the complex state vector in terms of the real coordi- 
nates. The operator J, that satisfies J 2 = —I, is called 
the complex structure. As it is apparent when comparing 
M-lUfy and ft4.1Sty , it implements the multiplication by the 
imaginary unity i in the real presentation of the Hilbert 
space. Observe that J and -Hr commute, this is due to 
the fact that the latter comes from an operator H in the 
complex Hilbert space. 

As the norm of a state can be written as the function 
associated to the identity operator, it is trivial to prove: 



Lemma 4 The dynamics defined by Eqs ^4.10\ l ( equiva- 
lently by preserves the norm of the state. 

Proof It is immediate to see that the evolution of the 
function /j is given by the Poisson bracket with the 
Hamiltonian: 

But because of the properties of the bracket 

{fH, A} = f[H,I] = 0. 

Thus the dynamics preserves the norm of the state, and 
the flow is restricted to the sphere Sq. □ 



V. EHRENFEST DYNAMICS AS A HAMILTONIAN 
SYSTEM 

In this section we show how to put together the dy- 
namics of a quantum and a classical system following the 
presentation of the two previous sections. We are describ- 
ing thus a physical system characterized by the following 
elements: 



8 



A. The set of states of our system 

• First, let H be a Hilbert space which describes the 
quantum degrees of freedom of our system. For ex- 
ample, it could describe the electronic subsystem; 
in this case, it is the vector space corresponding to 
the completely antisymmetric representation of the 
permutation group Sn (i.e., a set of Slater determi- 
nants), where N is the number of electrons of the 
system and each electron lives in a Hilbert space 
of dimension M. Thus, the dimension of H will be 

We know that it is a complex vector space, but we 
prefer to consider it as a real vector space with the 
double of degrees of freedom and denote it as Mq. 
Also, in correspondence with the Hilbert space vec- 
tors in the usual formalism of quantum mechanics, 
several states in Mq represent the same physical 
state. To consider true physical states one should 
extract only those corresponding to the projective 
space, which can be identified with a submanifold 
of Mq . A more general approach is to consider the 
sphere of states with norm equal to one, Sq, and 
take into account the phase transformations gener- 
ated by Eq (|4.ip in a proper way. We will discuss 
this in the following sections. 

• Second, let Mq be a differentiable manifold which 
contains the classical degrees of freedom. We shall 
assume it to be a phase space, and thus it will have 
an even number of degrees of freedom and it shall 
be endowed with a non degenerate Poisson bracket 
that in Darboux cordinates (R, P) reads 

r, , = df dg df dg 
U ' 9]C Y dPj dRj dRj dPj ■ 

• Third, we let our state space S be the Cartesian 
product of both manifolds, 

S = M c x M Q . 

Such a description has important implications: it 
is possible to consider each subsystem separately in 
a proper way but it is not possible to entangle the 
subsystems with one another. As long as Ehrenfest 
dynamics disregards this possibility, the choice of 
the Cartesian product is the most natural one. 

Example 2 If we consider a simple case, where we have 
one nucleus moving in a three dimensional domain and 
the electron state is considered to belong to a two-level 
system, the situation would be: 

• = (R,P,q 1 ,q2,pi,P2), (Qi,Q2,Pi,P2) e K 4 , 

where R represents the position of the nucleus, and P rep- 
resents its linear momentum. The tetrad (qi,q2,Pi,P2) 



represents the set of four real coordinates which corre- 
spond to the representation of the state of the two-level 
system on a real vector space ( of four 'real ' dimensions 
which corresponds to a 'complex' two-dimensional vector 
space). 

As a conclusion from the example above, we use as 
coordinates for our states: 

• The positions and momenta of the nuclei: 

(R,P)eM c . (5.1) 

We will have 3Nc + 3Nc of these, for Nq the num- 
ber of nuclei of the system. 

• The real and imaginary parts of the coordinates 
of the Hilbert space elements with respect to some 
basis: 

(q,p) E Mq. (5.2) 

We will have Nq + Nq of these, for Nq the complex 
dimension of the Hilbert space H. 

B. The observables 

To represent the physical magnitudes we must consider 
also the classical-quantum observables from a new per- 
spective. 

Our observables must be functions defined on the state 
space S = Mq x Mq. We know from our discussion in 
the case of a purely quantum system that any function 
of the form (|4.7p produces an evolution, via the Poisson 
bracket, which preserves the norm. In the MQCD case, 
we can easily write the analogue of the vector field (|4.1[) 
by writting: 

r Q = i®r. (5.3) 

This is again the infinitesimal generator of phase trans- 
formations for the quantum subsystem, but written at 
the level of the global state space Mc x Mq. A rea- 
sonable property to be asked to the functions chosen to 
represent our observables is to be constant under this 
transformation. From a mathematical point of view we 
can write such a condition as follows: 

Definition 5 We shall define the set of possible observ- 
ables, O, as the set of all C°° -functions on the set 
Mc x Mq which are constant under phase changes on 
the quantum degrees of freedom ,i.e., 

O = {f e C°°(M C xM q ) I r Q / = 0}. (5.4) 

As we shall see later, this choice reflects the fact that, 
when considered coupled together, the nonlinearity of 
classical mechanics expands also to MQCD. 

We would like to remark that because of the choice of 
the set of states as a Cartesian product of the classical 
states and the quantum states, we can consider as subsets 
of the set of observables: 
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• The set of classical functions: these are functions 
which depend only on the classical degrees of free- 
dom. Mathematically, they can be written as those 
functions / € O such that there exists a function 
fc e C°°(M C ) such that 

f(R,P,q,^) = fc(R,P). 

We denote this subset as Oc- An example of a 
function belonging to this set is the linear momen- 
tum of the nuclei. 

• the set of generalized quantum functions: 

functions which depend only on the quantum de- 
grees of freedom and which are constant under 
changes in the global phase. Mathematically, they 
can be written as those functions / € O such that 
there exists a function /q € C°°(Mq) such that 

f(R,P,q,p) = f Q (g,p); T(f Q ) = 0. (5.5) 

We denote these functions as Oq. We have added 
the adjective "generalized" because this set is too 
large to represent only the set of pure quantum 
observables. These later functions, should be con- 
sidered, when necessary, as a smaller subset, which 
corresponds to the set of functions defined in Eq. 
(|4.7|) . We denote this smaller subset as Oq. An 
example of a function belonging to Oq is the linear 
momentum of the electrons. 

• A third interesting subset is the set of arbitrary 
linear combinations of the subsets above, i.e., those 
functions which are written as the sum of a purely 
classical function and a purely quantum one: 

f(R,P,q,p) =f C (R,P) + f Q (q,p). (5.6) 

We will denote this set as Oc+q- An element of 
this set of functions is the total linear momentum 
of the composed system. 

We would like to make a final but very important re- 
mark. We have not chosen the set of observables as 

V = [f A e C°°(M C x Mq)\Sa = (ip\A(R,P)ip)} , 

(5.7) 

for A(R, P) a linear operator on the Hilbert space T-L 
depending on the classical degrees of freedom because of 
two reasons: 

• It is evident that the set above is a subset of (|5.4[) 
and thus we are not loosing any of these observ- 
ables. But it is a well known property that Ehren- 
fest dynamics is not linear and then if we consider 
the operator describing the evolution of the system, 
it can not belong to the set above. We must thus 
enlarge the set (|5.7j) . 

• We are going to introduce in the next section a 
Poisson bracket on the space of observables. For 
that bracket to close a Poisson algebra, we need to 
consider the whole set (|5.4[) . 



It is important to notice that in the set (|5.4[) there 
are operators which are not representing linear operators 
for the quantum part of the system and hence the set of 
properties listed in Section HVl for the pure quantum case 
are meaningless for them. But this is a natural feature of 
the dynamics we are considering, because of its nonlinear 
nature. 



C. Geometry and the Poisson bracket on the 
classical-quantum world 

Finally, we must combine the quantum and the clas- 
sical description in order to provide a unified descrip- 
tion of our system of interest. As we assume that both 
the classical and the quantum subsystems are endowed 
with Poisson brackets, we face the same problem we have 
when combining, from a classical mechanics perspective, 
two classical systems. Therefore it is immediate to con- 
clude that the corresponding Poisson structures can be 
combined as: 

{;-} = {;-}c + h- 1 {;-} Q , (5.8) 

where the term {•, -}c acts on the degrees of freedom 
of the first manifold and {•, - }q acts on the degrees of 
freedom of the second one. It is a known fact^ that such 
a superposition of Poisson brackets always produces a 
well defined Poisson structure in the product space. 

The set of pure classical functions Oc and the set of 
quantum generalized functions Oq are closed under the 
Poisson bracket. The same happens with the quantum 
functions Oq and the set of linear combinations Oc+q ■ 
In mathematical terms, what we have is a family of Pois- 
son subalgcbras. This property ensures that the descrip- 
tion of purely classical or purely quantum systems, or 
even both systems at once but uncoupled to each other, 
can be done within the formalism. 

Once the Poisson bracket on Mc x Mq has been intro- 
duced we can express again the constraint we introduced 
in the definition of the observables in Poisson terms. 
Thus we find that in a completely analogous way to the 
pure quantum case, we can prove that 

Lemma 5 The condition in Eq. (|5.4[) 

r<j(/) = o 

is equivalent to ask the function f to Poisson- commute 
with the function fj — J^kilk + Pfc)> *- e -> 

r Q (/) = o^{/i,/} = o. 



D. The definition of the dynamics 

From Sections IIIII and IIV1 our formulation of MQCD 
can be implemented on: 
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• The manifold which represents the set of states by 
defining a vector field whose integral curves rep- 
resent the solutions of the dynamics (Schrddingcr 
picture). 



• The set of functions (please note the differences 
between the classical and the quantum cases) de- 
fined on the set of states which represent the set 
of observables of the system. In this case the Pois- 
son bracket of the functions with the Hamiltonian 
of the system defines the corresponding evolution 
(Heisenberg picture). 



Both approaches are not disconnected, since they can 
be easily related: 



Xh — {/h, ■}, 



(5.9) 



where wc denote by Xh the vector field which represents 
the dynamics on the phase space and by fu the function 
which corresponds to the Hamiltonian of the complete 
system. 

We can now proceed to our first goal: to provide a 
Hamiltonian description of Ehrcnfcst dynamics in terms 
of a Poisson structure. We thus define the following 
Hamiltonian system: 



• A state space corresponding to the Cartesian prod- 
uct Mc x Mq. 



As a result, the dynamics of both subsystems are ob- 
tained easily. In the Schrodingcr picture we obtain: 



R 
P 



dfn = M _ ip 



= h 



dP 

dfn = 
dR ' 
idf H 



Pi = -ft 



dpi ' 
-idf H 
dqi ' 



PNr 



n 
-ft 



i df H 



9pn q ' 

-1 dfn 



dqN Q 



(5.11) 

v a d((^(q,p)\H e (R)^,p))), (5.12) 

(5.13) 
(5.14) 

(5.15) 
(5.16) 



Remark that, form a formal point of view, this is anal- 
ogous to Example [TJ This set of equations corresponds 
exactly with Ehrenfest dynamics. 

The final point is to prove the following lemma: 

Lemma 6 The dynamics preserves the set of observables 
O. 

Proof From Lemma |4j an observable belongs to O if it 
Poisson-commutes with f\. Thus, as fu £ O, if we con- 
sider an observable / G 0, by the Jacobi identity: 

{h, Uh, /}} = -{/, {ft, .M} - {f H , {/, h}} = o. 

(5.17) 
□ 



VI. STATISTICS OF THE EHRENFEST DYNAMICS 



• A set of operators corresponding to the set of func- 
tions O defined in Eq. (|5.4j) . 



• The Poisson bracket defined in Eq. (|5.8|) . 



• And finally, the dynamics introduced by the follow- 
ing Hamiltonian function: 

f H (R,P,q,p) = + (i>(q,p)We(RMq,P)}, 

(5.10) 

where H e is the expression of the electronic Hamil- 
tonian, Mj are the masses of the classical subsys- 
tem of the nuclei and ?p(q,p) is the real-space rep- 
resentation of the state ip analogous to Eq. (|4.14l) . 



The next step in our work is the definition of a sta- 
tistical system associated to the dynamics we introduced 
above. The first ingredient for that is the definition of 
the distributions we shall be describing the system with. 

The main conclusion from the previous section is that 
Ehrenfest dynamics can be described as a Hamiltonian 
system on a Poisson manifold. Therefore, we are in a 
situation similar to a standard classical system. We have 
seen that the dynamics preserves the submanifold Mc x 
Sq (Sq being the sphere in Eq. (|4.2p ). and thus it is 
natural to consider such a manifold as the space of our 
statistical system. 

We have thus to construct now a statistical system 
composed of two subsystems. We may think in introduc- 
ing a total distribution factorizing as the product of a 
classical and a quantum distribution. But as the events 
of the classical and the quantum regimes are not inde- 
pendent, the two probabilities must be combined and not 
defined through two factorizing functions. 

Hence we shall consider a density Fqc £ O and the 
canonical volume element dy,Qc in the quantum-classical 
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phase space, restricted to Mq x Sq: 

d^Qc = dficdClQ, (6.1) 

where d[ic is given by (|6.3[) and d£lQ is the volume ele- 
ment on a (2Nq — l)-dimcnsional sphere (see below). 
With these ingredients we introduce the following 

Definition 6 Given an observable M £ O, we define the 
macroscopic average value, (M), as 

(M)= [ M(R,P,q,p)F Q c(R,P,q,p)d»QC (6.2) 

J M c xSq 

In the rest of the section we shall discuss in more detail 
the ingredients of this definition. 

We start by the volume form. Given the classical phase 
space Mc with Darboux coordinates (R, P) (they always 
exist locally 4 ^) we define the volume element in the clas- 
sical phase space by 

d f i c =Y[dPjdR j7 (6.3) 
J 

that, as it is well known, is invariant under any purely 
classical Hamiltonian evolution. One may wonder to 
what extent the volume form depends on the choice of 
coordinates. But taking into account that two different 
sets of Darboux coordinates differ by a canonical trans- 
formation, and using the fact that these transformations 
always have Jacobian equal to one, one concludes that 
the volume form does not depend on the choice of coor- 
dinates. 

We might proceed in the same way for the quantum 
part of our system, as we have written the quantum dy- 
namics like a classical Hamiltonian one and, therefore, 
we could use the canonical invariant volume element de- 
scribe in the previous paragraph dfiQ on Mq. However, 
new complications appear in this case as in the quantum 
part of the phase space we restrict the integration to Sq 
and then we look for an invariant volume form in the unit 
sphere, rather than in the full Hilbcrt space. 

The new ingredient that makes the restriction possible 
is the fact that all our observables and of course, also the 
Hamiltonian, are killed by T. This implies, as discussed 
before, that f\ is a constant of motion. Using this we 
decompose 

d^Q = dfjdh Q , 

which is possible on % \ {0}. The decomposition is not 
unique, but the restriction (pullback) of dQ,Q to the unit 
sphere gives a uniquely defined, invariant volume form 
on Sq that we will denote by df2g. The decomposition 
is nothing but the factorization into the radial part and 
the solid angle volume element. 

Example 3 For the simple case we studied above, where 
% = C 2 and hence Mq = K 4 and Sq is a three dimen- 
sional sphere, the volume element above would be: 

dfl Q = qidq 2 dpidp 2 - q 2 dqidpidp 2 

+pidqidq 2 dp 2 - P2dqidq 2 dpi. (6.4) 



Now, we finally put together the two ingredients to 
obtain the invariant volume element dfiQc in cq. (j6.ll) . 

We discuss now the properties we must require from 
Fqc in order for Eq. (|6.2p to correcly define the statisti- 
cal mechanics for the Ehrenfest dynamics. 

• The expected value of the constant observable 
should be that constant, which implies that the in- 
tegral on the whole set of states is equal to one: 

f F QC dfi QC = 1. (6.5) 

Jm c xSq 

• The average, for any observable Ja of the form (|5.7[) 
associated to a positive definite hermitian operator 
A should be positive. This implies the usual re- 
quirements of positive probability density in purely 
classical statistical mechanics. 

So far we have discussed general considerations about 
the statistical mechanics of our system. Now wc should 
include the dynamics in the game. Given that we have 
formulated the Ehrenfest dynamics as a Hamiltonian sys- 
tem, it is well known^ which is the evolution of the den- 
sity function Fqc- It is given by the Liouville equation 

^ + {f H ,F QC } = 0, (6.6) 

where, in the derivation of the equation, it is an essential 
requirement the invariance of the volume form under the 
evolution of the system. 

The equilibrium statistical mechanics is obtained by 
requiring Fqc = 0, which making use of (|6.6[) is equiva- 
lent to 

{f H ,F QC } = 0, 

or, in other words, Fqc should be a constant of motion. 
An obvious non trivial constant of motion is given by any 
function of fn- The question is if there are more. 

To answer this we should examine closely our particu- 
lar dynamics. In several occasions we have mentioned the 
non linear character of a generic Ehrenfest system. This 
is an essential ingredient, because due to this fact we do 
not expect any other constant of motion that a function 
of the Hamiltonian itself. This is in contrast with the 
case of linear equations of motion, where the system is 
necessarily integrable. 

The question of the existence of more constants of mo- 
tion and, therefore, more candidates for the equilibrium 
density Fqc can be elucidated with a simple example 
that we discuss below. In it we show the presence of 
crgodic regions (open regions in the leave of constant en- 
ergy densely covered by a single orbit). This rules out 
the existence of constants of motion different from the 
constant function in the region (i.e., functions of fa)- 

Example 4 In the following example we will study a 
simple toy model in which the coupling between classi- 
cal and quantum degrees of freedom gives rise to chaotic 
behavior and the appearence of ergodic regions. 
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FIG. 1. The plot shows a single orbit of the Poincare map of 
Eq. (|6.7|) at 8 = 0. The angle (/> £ [—ir,ir] is plotted in the 
horizontal axis and I& £ [0, 1] in the vertical one. We have 
taken e = 0.8. The energy determines Ig and its actual value 
does not affect the dynamics of <j> and Is- 



The system consists of a complex two dimensional 
Hilbert space Mq = C 2 and a classical 2-D phase space 
where we define a 1-D harmonic oscillator. Using co- 
ordinates (le,0) for the classical variables (action- angle 
coordinates for the oscillator) and f G C 2 me define the 
following Hamiltonian 

fa =Ie + + ecos(0)<r a! |tf), 

with o~ x , o~ z the Pauli sigma matrices. 

We parametrize the normalized quantum state by 



e«Jl-I$ 



One can now easily check that 1$ and <f> are canonical 
conjugate variables. 

In these variables the Hamiltonian reads 



f H =h 



el^/l- /|cos(0) cos(</>), 



(6.7) 



where e measures the coupling of the classical and quan- 
tum systems. 

In the limit of vanishing e the system is integrable and 
actually linear in these coordinates. However for non 
vanishing e the model becomes non linear and, as we will 
show below, regions of chaotic motion emerge. 

In order to understand the behaviour of our system 
it rs useful to study Us Potncare map (see Ref. ft for 
the definition). To this end we take the transversal (or 
Poincare) section at 8 = and, taking into account the 
conservation of energy, only two coordinates are needed 
to describe the map, we have chosen the quantum vari- 
ables 1$ and <f>. 

By examining the plot of Fig. [7J one sees that the orbit 
densely fills a large region of the Poincare section and 
therefore any constant of motion should be the constant 
function in that region. In the presence of such an ergodic 
evolution the only constants of motion are functions of 
the Hamiltonian. 



As a consequence we can claim that, genericaly within 
the Ehrenfest dynamics, the only functions which com- 
mute with the Hamiltonian in (|6.6|) arc those which are 
function of the Hamiltonian fjj itself. As the Poisson 
bracket is skew-symmetric, the property follows trivially. 
Ergodicity helps us to assert, roughly speaking, that in 
general no other function in O will be a constant of the 
motion. 

Thus, at this stage, we can finish our construction with 
the equilibrium distribution associated to Ehrenfest dy- 
namics. Taking into account the ergodicity of the dy- 
namics, we can impose the equal-probabilities condition 
to the configurations of our isolated system, which leads 
to the microcanonical ensemble, (see Ref. |43I ). Then, the 
probability density for the canonical ensemble turns out 
to be: 



FS c (R,P,q,p) = Z 



(6.8) 



where Z = J AIcxSq s "^ H d/j,Qc, as it can be seen, in 
a standard derivation, for example in Refs. l36l . l37l . and 
[iil . where being proportional to the temperature 

of the isolated system, is the parameter governing the 
equilibrium between the different parts of the system. An 
information-theoretic approach to the equiprobability in 
the microcanonical ensemble and to Eq. (|6.8p may be 
seen in Ref. or in Refs. |46| and H3 for a much more 
modern and mathematically sound exposition. 



VII. CONCLUSIONS AND FUTURE WORK 

In this paper we have constructed a rigorous Hamilto- 
nian description of the Ehrenfest dynamics of an isolated 
system by combining the Poisson brackets formulation 
of classical mechanics with the geometric formulation of 
quantum mechanics. We have also constructed the cor- 
responding statistical description and obtained the asso- 
ciated Liouville equation. Finally, after verifying numer- 
ically that the Ehrenfest dynamics is ergodic, we justify 
the equilibrium distribution produced by it. 

The definition of the Hamiltonian in Eq. (|5.10[) and 
its associated canonical equilibrium distribution (|6.8[) al- 
lows us now to use (classical) Monte Carlo methods for 
computing canonical equilibrium averages in our system, 
given by the expression (|6.2D . 

Also it is straightforward to extend the non-stochastic 
method proposed by Nos 6 16 i 17 to our Hamiltonian (|5.10|) . 

However, if we wish to perform MD simulations us- 
ing stochastic methods, we should first construct the 
Langcvin dynamics associated with our Hamiltonian 
(|5.10[) and equations (|5.11|) - (|5.16[) . To develop this pro- 
gram, we also need an analogue of the Fokkcr-Planck 
equation associated with those Langevin equations and 
then we have to check whether its solution at infinite 
time approaches (|6.8[) . Presently we are working on this 
points^, i.e., on the extension of our formalism to an as- 
sociated stochastic MD that can be used to rigorously 
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probe non-equilibrium phenomena at constant tempera- 
ture. 
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